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ABSTRACT 



We investigate statistical equilibrium of neutral and singly-ionized strontium in late-type stellar atmospheres. Particular attention is 
given to the completeness of the model atom, which includes new energy levels, transition probabilities, photoionization and electron- 
impact excitation cross-sections computed with the R-matrix method. The NUIE model is applied to the analysis of Sr I and Sr II 
lines in the spectra of the Sun, Procyon, Arcturus, and HD 122563, showing a significant improvement in the ionization balance 
compared to LIE line formation calculations, which predict abundance discrepancies of up to 0.5 dex. The solar Sr abundance is 
log e = 2.93 ± 0.04 dex, in agreement with the meteorites. A grid of MUTE abundance corrections for Sr I and Sr II lines covering a 
large range of stellar parameters is presented. 

Key words, atomic data - line: formation - radiative transfer - sun: abundances - stars: abundances 



1. Introduction 

Spectroscopic observations of low-mass stars have shaped our 
understanding of Galactic evolution and stellar nucleosynthe- 
sis. Strontium, as one of the abundant r-process elements, has 
been extensively investigated in the past few decades. However, 
its main production site has not yet been identified: the ob- 
served abundances of Sr in metal-poor stars are far too large 
to be explained by conventional rapid neutron-capture nucle- 
osynthesis in SNe II, suggesting some alternativ e exotic scenar- 
ios, s uch as the light element primary process (iTravaglio et al.l 
|2004|) , rp-process in acc retion disks around low-mass black 
holes ("Schatz et al.' 1998!), black hole - neutron star mergers 



dSurman et al. 2008). high-entropy winds in SN II ( Farouai et al 
20101) . and low-mass electro n-capture supernovae dWanajo et al 



201 ih to name just a few (see lJose &Iliadisl20lil and references 
therein). 

Until recently, determinations of Sr abundances in metal- 
poor stars relied almost exclusively on the two near-UV lines 
of Sr II, which are sufficiently strong to be detected also in the 
spectra of moderate-to-low resolution and signal-to-noise. The 
drawback is that in spectra of stars typically used for studies of 
Galactic chemical evolution, [Fe/H] > -2, these lines saturate 
and develop pronounced damping wings, overlapping with var- 
ious atomic and molecular blends. Thus, at higher metallicities 
a preference is sometimes given to the weak Sri line at 4607.34 
A and/or the Sr ii line at 4 1 6 1 A. However, there is evidence that 
the Sr I lines may be subject to non -local thermodynamic equi - 
librium (hereafter, NLTE) effects dBarklem & O'Maral I2OO0I) . 
which has been supported by ab initio calculations solving for 
radiative tr ansfer in NLTE for a small grid or red giant model at- 
mospheres dShort & Hauschild3l2006l) . In a few studies utilizing 
near-IR spectra, als o the Srn triplet (10037 , 10 327, and 10915 
A have been used dAndrievskv et al.ll201 Ih . In the majority of 
published studies, the preference is given to one ionization stage 



only (e.g.. lJehin et al.lll999|). and only a few studies investigated 
both Sri and Srii hnes dGratton & Snedenlll994t ICowan et aP 
l2002h . finding discrepant results. 

In this study, we perform for the first time a NLTE analy- 
sis of the Sri and Srii lines in spectra of late-type stars. The 
new atomic model was constructed from the state-of-art atomic 
data, computed specifically for this work. The NLTE model atom 
is tested on a number of reference stars with parameters deter- 
mined by other independent methods. Furthermore, we present a 
large grid of NLTE abundance corrections for Sri and Srii lines. 
The results presented in this work will be applied to the anal- 
ysis of a representative sample of metal-poor stars observed at 
a very-high resolution and signal-to-noise in Hansen et al. (in 
prep.) Furthermore, we plan to undertake the NLTE Sr abun- 
dance analysis o f the thick-disk and h alo stars with spectroscopic 
parameters from lRuchti et alj d201 lb . 

2. Methods 

2.1. Model atmospheres and NLTE calculations 

All calculations in this work were performed with classical ID 
LTE plane-parallel model atmospheres. We used the MAFAGS- 
OS (Grupp 2004abi) and MAFAGS-ODF models, which are 
well-adapted for the analysis of late-type stars. A brief descrip- 
tion of these models a nd comparison with ot her models of a 
similar type (MARCS, iGustafsson et al.ll2008h is presented in 
iBergemann et al.l d2012h . 

The NLTE statistical equilibrium calculations were per- 
formed with the revis ed version of the DETAIL code 
dButler & Giddingslll985h . The statistical equilibrium and radia- 
tive transfer equations are so lved by the accelerated A-iteratio n 
method in the formulation of lRvbicki & Hummed dl99lLll992h . 
The method allows for self-consistent treatment of overlapping 
lines and continua. A description of the code with some recent 
modifications related to the treatment of background line opacity 
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can be found in lBergemann et aP (l2012h . In the statistical equi- 
librium calculations, we treat the diagnostic lines investigated 
here with Voigt profiles, whereas all other Sr lines were com- 
puted with a Gaussian profile with 13 frequency points. 

The LTE and NLTE abundances for the reference stars, as 
well as the NLTE abundance coiTections, were determined by 
full s pectrum synt hesis with the revised version of the SIU 
code (lReetzlll999h . which was adapted for automated NLTE 
abundance calculations (L. Sbordone, private communication). 
The line lists in SIU have bee n continuo usly updated by the 
members of the Munich group dGehren et a l. 2004, 200ct and 
they are mainly based on the Kurucfl Hannover 0, and NIST 
(iRalchenko et alJl2012h compilations. 

2.2. Model atom 

The model atom of Sr was constructed using the available 
atomic data from the Hannover and NIST databases and sup- 
plemented by our new calculation^ of atomic energy levels, 
dipole-permitted and forbidden transitions (see below). Thus, 
the model is the most complete representation of atomic sys- 
tem of Sri and Srn, to date. In particular, we also include 135 
Srn transitions fr om the Cou lomb approximation calculations 
of ILindgard & N ielsonI (Il977h . The neutral atom is, thus, repre- 
sented by 141 levels with the principle quantum numbers n* < 
20, and the uppermost level 20f 'F° is located at 5.66 eV, i.e. 0.03 
eV from the 1-st ionization limit at 5.69 eV. The singly-ionized 
atom includes 49 levels including the highest experimentally- 
observed level 14d ^D, which is located at 10.684 eV. The model 
is closed by the Sr iii ground state 4p^ ' So . The number of dipole- 
allowed transitions is 336 and 214 for Sri and Srii, respectively. 



Our model of Sr is similar to that of lAndrievskv et aP (1201 ih 
with respect to the term structure and the number of dipole- 
permitted transitions of Srii. However, we include more ac- 
curate transition probabilities for Srii. More importantly, our 
model atom includes new quantum-mechanical photoionization 
cross-sections for Sri and e" impact excitation cross-sections for 
dipole-permitted and forbidden transitions in Srn ( s ee bel ow). 
The model atoms employed b y IShort & Hauschild i (l2006l) and 
iMashonkina & GehrenI (1200 ih are much simpler than that em- 
ployed in this work (84 and 41 total levels, respectively). 

2.2.1 . New calculations of atomic data 

We compute single-electron orbitals for the t arget ion s ^ Sr^ 
with the atomic structure code autostructure by iBadnelll ( 19861 
Il997h. This code is an extension of the program superstruc- 
ture dEissner & Nussbaumei1ll969h . computes fine-structure CI 
energy levels, and radiative and Auger rates in a Breit-Pauli 
relativistic framework. Single electron orbitals are constructed 



' kurucz.harvard.edu/ 

2 www.pmp.uni-hannover.de/cgi-bin/ssi/test/kurucz/s ekur.htmI. The 



atomic line data used in this database are taken from iKurucz & BeiH 
( fl995l) . 

^ The new atomic data for Sr presented in the paper and in Bautista 
et al. (2002), including the energy levels, /-values, photoionization and 
e-impact excitation cross-sections, can be provided upon request. 

* We note that different notations are standard for atomic physics 
and spectroscopy, and we chose to respect this difference in our work. 
Ionization stage of an element is thus indicated by a numerical super- 
script in Sect. 2.2.1 describing the atomic physics calculations and by a 
roman numerical in astrophysical application in the further sections. 





'S ¥ V 'd 'd° V 'S 'P° 'D 'D° 'F° 



Fig. 1. Grotrian diagram of the Sr i (bottom) and Sr ii (top) model 
atoms. Only dipole-allowed transitions are shown. This is the 
most complete model atom utilized for NLTE calculations of 
FGK stars so far. 



by diagonalizing the non-relativistic Hamiltonian, within a sta- 
tistical T homas-Fermi-Dirac-Amaldi (TFDA) model potential 
(lEissner & Nussbaumer 1969). The scaling parameters are op- 
timized by minimizing a weighted sum of the LS term energies. 
We employ a very extensive configuration expansion with con- 
figurations of the form Ap^nl, with n going from 4 to 6 and / < 3, 
and configurations with multiple promotion from the 4p orbital, 
i.e. ApiSs'M'nl with 3 < ^ < 6, < r < 2, and < f < 2. These 
multiple promotions from the 4p orbital are essential in comput- 
ing the photoionization cross sections and that will be the subject 
of a future publication. 

Our target representations give term energies in reasonable 
agreement, within ~ 5 percent with experimental values. 

The photoio nization cross-sect ions are computed with the R- 
matrix method (iBurke et al.lll97l|) . The present calculation for 
Sr" includes the lowest 67 LS-terms of the target in the close 
coupling expansion and all short range (N-i-l)-electron configu- 
rations that results from adding an electron to the target config- 
urations. We find 197 singlet and triplet bound terms of Sr" with 
n < 9 and L < 5, and compute photoionization cross sections 
for all of them. The cross sections are computed at an even en- 
ergy mesh of 5 X 10"'' Ry from threshold up to 0.23 Ry, mesh 
of 10"^ Ry from 0.23 Ry to 2.3 Ry (roughly the highest tar- 
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Fig. 2. Partial photoionization cross-sections for the Sr i ground Fig. 3. Total photoionization cross-sections for the selected Sr i 
state into the lowest six states of the Srii target. The cross- levels. Dotted lines show hydrogenic cross-sections. The new 



sections are dominated by prominent resonances. 



cross-sections are up to few orders of magnitude larger than the 
commonly-used hydrogenic approximation. 



get threshold), and another mesh of 0.02 Ry from 2.3 Ry to 4.3 
Ry. Figure |2] shows the partial cross sections for ionization from 
the ground term of Sr" into each of the first six terms of the 
target Sr^. We note, however, that DETAIL is not yet able to 
consistently treat ionization to specific states of the target ion. 
Therefore, total photoionization cross-sections are adopted here. 
Figure [3] shows the total cross sections for the lowest six excited 
state of Sr". For comparison, we show the hydrogenic cross- 
sections computed with effective principal quantum number The 
latter reproduce the functional dependence of the cross-section 
with energy. However, the quantum-mechanical photoionization 
cross-sections are typically larger both in the background (fac- 
tor of two to five) and in the resonances (up to three orders of 
magnitude). 

2.2.2. Collision rates 

Electron impact excitation rate coefficients for the 49 levels of 
Sr II levels were taken from lBautista et all (l2002h . The data were 
computed with the same technique, i.e., the close-coupling ap- 
proximation with the R-matrix method. For all other levels of 
Sri and Srii th e electron collision rate s were computed using 
the formulae of Ivan RegemorteJ (1 19621) for dipole-allowed and 
lAllenl(ll973h for dipole-forbidden transitions. Cross-sections for 
transitions caused by inelastic collisions with Hi atoms are 
basically unknown for any atom heavier than Mg. The only 
available formula developed originally for collisions between 



equal H-like atoms ("Drawin" 1 968', ' 1 969t) and later slightly mod- 
ified by Steenbock & Holweger ( 1984} and lLamber t ( 1993i) was 
shown to over-estimate the rates of bound-bound transitions by 
two to seven orders of magnitude (IBelvaev & Barkle"ml 120031 : 
iBarklem et a02012l) . Furthermore, the charge transfer processes 
can not be described by this simple classical formula at all. Our 
tests with various scaling factors to the Drawin's Hi inelastic 
collision cross-sections demonstrated that excitation and ioniza- 
tion balance of Sri/Srii in the reference stars can be satisfied 
simultaneously only if the efficiency of Hi collisions is very 
low (Sect. 4), few orders of magnitude lower than that given 
by the Drawin's recipe, which is consistent with the quantum- 
mechanical results mentioned above. On these grounds we do 
not include inelastic H i collisions in our reference model atom. 

2.3. Line selection and atomic data 

In spectra of FGK stars, only a few Sr lines are useful for abun- 
dance determinations. These are two Sr i and 6 Sr ii lines. The Sr i 
resonance line at 4607 A is weak, and, thus, sufficiently reliable 
for more metal-rich stars ([Fe/H]> -1.5) observed at high reso- 
lution. The two resonance Sr ii lines, in contrast, remain strong 
even in spectra of very metal-poor stars, however, severe blends 
in the inner and outer wings lead to systematic errors in abun- 
dances, when the blends are not properly accounted for The 
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Table 1. Lines of Sri and Srii selected for abundance calcula- 
tions. 



A 
A 



[eV] 



low 



'-up 

[eV] 



up logg/ log(r/A'H)'^ logCs^ 

s" 



Sri 

4607.33 
7070.01 
Srii 

4077.71" 

4167.80 

4215.52° 

10036.65 

10327.31 

10914.89 



0.00 5s- 'So 
1.85 5p3p° 



2.69 
3.60 



5p'Pi 0.283 
6s 3S, -0.020 



0.00 5s ^S 

2 



1/2 

2.94 5pX,2 
0.00 Ss^Si/j 
1.81 4d2D,„ 



1.84 4d-D^ 
1.80 4d2D ' 



2po 

5.92 6s2Si/2 

2po 

>o 

3/2 
3/2 
1/2 



3.04 5p .3/2 



2.94 5p . 
3.04 5p2p° 
3.04 5p2p; 
2.94 5p2p' 



0.158 
-0.502 
-0.155 
-1.194 
-0.240 
-0.474 



-7.53 
-7.15 

-7.81 
-7.81 
-7.81 
-7.64 
-7.64 
-7.64 



-31.2 
-30.2 

-32.0 
-32.0 
-32.0 
-31.5 
-31.5 
-31.5 



Notes. 



blended. <*' Log of the FWHM per H atom at 10000 K. '-'^ Ce are in 



the units of cm s 



Table 2. Hyperfine structure and isotopic shift for the 4077 A 
line. 



Isotope 


A 
A 




87 


4077.697 


-1.645 


87 


4077.699 


-1.485 


84 


4077.708 


-2.094 


86 


4077.709 


-0.848 


88 


4077.710 


0.075 


87 


4077.724 


-1.465 


87 


4077.725 


-1.956 



three near-IR Srn lines at Ifim appear to be un-blended and are 
sufficiently strong to be detected even at [Fe/H] < -2. 

The lines selected in this work are given in Table [T] The 
transition probabilities were taken from different experimental 
and theoretical sources. According to NIST, the Sri ^/-values 
are very accurate, the uncertainty is lower than 1 percent for the 
4607 A line and better than 10 percent for the 7070 A subordi- 
na te line. These g f- valu es are taken from the labor atory analysis 
of lParkinsonI ( Il976l) and lGarcia & CamposI ( Il988h . respectively. 
The oscillator strengths of the near-IR Srn lines were adopted 
from the only available laboratory results of Gallagher (1967). 
These are fully consistent with the most recent ab initio calcu- 
lations of iBautista et al.l (l2002h . which are also adopted in the 
NLTE model atom. 

Damping widths for the calculation of broadening due to 
elastic collisions with Hi are available f or the Srn lines from 
the quantum-mechanical calculations of iBarklem et alj (|2000|) 
and for the Sri lines they were kindly provided by the ref- 
eree. We adopt the values from Barklem et al. for the near-IR 
Srii transitions. However, for the strong resonance Srii lines, 
which are very sensitive to this parameter, our spectrum syn- 
thesis calculations indicate somewhat lower values, by ~ 20 per- 
cent. The question is whether this difference can be attributed 
to the uncertainties of the theoretical data arising because of 
certain approximations in the Anstee, Barklem, and O'Mara 
(ABO) theory. These include representation of the interaction 
potential and collisional dynamics and have been investigated 



by iKerkeni et al.l (l2004l) for neutral atoms of Mg, Sr, Ca, and 
Na. They find that whereas the semi-classical description of a 
collision is a sufficiently-good approximation, the ABO poten- 
tials become rather inaccurate at small interatomic separations. 
However, the effect of the latter is rather to under-estimate the 
line width that clearly does not explain our finding. For the Srii 
4167 line we assume the same y as for the resonance lines, but 
note that the line is too weak even in the solar spectrum to be 
sensitive to this parameter at all. 

The adopted values of damping widths at 10000 K per H 
atom, y/Nu, as well as the commonly-used interaction constants 
in the van der Waals-type interaction Cg computed for T = 5780 
K, are given in th e Table [T] Here, the latter is a parameter in 
the lUnsoldl d 19551 Eq. 82,48) formula needed to reproduce the 
correct line width y at the given temperatur^. 

For the 4077 Sr n line, which will be used for the abundance 
determinations in Sect. 4, we included the hyperfine structure 
and isotopic shifts. The m agnetic dipole and elect ric quadrupole 
constants were taken from lBuchinger et"an (Il990l) : A (5s^ 'Sq)= 
-1000.5 MHz, A {5p-P°y,) = -36 MHz, B (5p2p°^2) = 88.5 



3/2' 



M Hz. The data are in agreement with the older experiment 
by iBorghs et alj ( 1 1983b and theoretical calculations by Yu et 
al. (2004). Isotope shifts for ^"^Sr, '^''Sr, '^^Sr were taken from 
iBorghs et alj (Il983h and the solar isotopic ratios were adoptecfl 
The wavelengths and logg/ values of the hyperfine structure 
(HFS) components are given in Table |2l We note that the 4215 
A Sr I line is only used for a comparative LTE to NLTE analysis 
for the reasons discussed in Sect. 3.3. 

3. Statistical equilibrium of Sr 

The departure coefficient^ for selected levels of Sr are shown 
as a function of continuum optical depth logT^ at 500 in Fig. 
|4] for the model atmospheres with parameters corresponding to 
dwarfs and giants with [Fe/H] = and -2.4. Only the levels, 
which give rise to the spectral lines selected for the abundance 
analysis (Table[T]i, are included in the plots. These are two levels 
of Sr I, forming the resonance lines at 4607 A, and two levels of 
Sr n, which are connected by the 4077 A transition commonly 
used in the analysis of metal-poor stars. 

3.1. NLTE effects in Sr\ 

Sr I with the ionization potential 5.69 eV is a trace atom in the the 
atmospheres of late-type stars. The ratio of Sri/Srn falls below 
10"^ above continuum optical depth unity, and the departures 
from LTE in the distribution of atomic level populations are al- 
most entirely due to overionization. The behavior of the Sr i de- 
parture coefficients with stellar paramet ers, thus, resembles tha t 
of the s imilar trace atoms, such as Co i dBergemann et alj|2010h 
or Cri dBergemann & Cescuttill20l"ol) . The levels are underpop- 
ulated and their /7/-factors monotonously decrease outwards. 
The ground state 5s^ 'S with the ionization edge at 2177 A is 
also strongly over-ionized and it decouples from the other low- 
excited levels because collisions can not bridge the large energy 



^ Note that Cg is often computed using the Unsold (1955, Eq. 82,55) 
approximation, in which the interaction constant derives from the effec- 
tive principal quantum number of lower and upper levels of a transition. 
However, t his approximat ion is known to yield too small damping con- 
stants (e.g. lBarkleml l200'^. and references therein) 

* http://www.nist.gov/pml/data/comp.cfm 

' The departure coefficients are defined as the ratio of NLTE to LTE 
level number densities, bj = ;j^™/;i™ 
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Fig. 4. Departure coefficients bi of selected Sri and Srii levels for different stellar parameters, T^ff, \ogg, and [Fe/H] (indicated in 
each sub-plot). Red: Sri levels, Ss^ 'S (soUd), 5p^P° (dotted), 5p 'P° (dashed). Black: Srii levels, 5s ^Sj^j (sohd), M^G^^^ (dotted), 
5p2p°^2 (dashed). 



gap of ~ 2 eV between them. Line transitions have a small influ- 
ence on the statistical equilibrium of Sri causing some changes 
of the level populations only in the solar-metallicity models. For 
example, in the metal-rich models (Fig.|4] left panel), the depar- 
ture coefficient of the 5p 'P° level drops below that of the 5s^ 
'S at logTc ~ -1, marking the depth at which the 4607 reso- 
nance line becomes optically thin, and spontaneous transitions 
depopulate the upper level of the transition. 

In the metal-poor models (Fig. H) right panel), the coupling 
of the levels by radiative bound-bound transitions and collisions 
diminished even though the collision rates between the low- 
excited levels are large enough to keep them are in detailed bal- 
ance even in the models with [Fe/H] = -1. The populations are 
now set primarily by the balance between radiative ionizations 
and recombinations, and the run of i?,-factors with log sim- 
ply reflects the mean intensity vs Planck function inequalities, 
Jv + By, at the frequencies of the level photoionization continua 
at each depth point. For example, in the model with Teff = 6000 
K, log^ = 4.2, and [Fe/H] - -2.4, the lowest energy levels 5p 
^P° and 5p 'P° decouple from each other and from the contin- 
uum at log Tc ~ 0. Then, in a very narrow interval of optical 
depths, -1 < log Tc < 0, their departure coefficients demonstrate 
a sudden drop following the steep temperature gradient of the at- 

* Note that we do not include the H i collisions as discussed in Sect. 

4. 



mosphere, and thus, rapidly increasing - By imbalance (note 
that matter becomes transparent to the radiation above ioniza- 
tion edges of both levels, 3202 and 4126 A at logr^ = +0.11 
respectively log Tc = +0.05). Further out, log Tc < -1, where the 
mean intensity as well as the local kinetic temperature remain 
roughly constant, the departure coefficients smooth out, slowly 
decreasing outwards. 

The NLTE effects on the formation of the Sr i line at 4607 
A are generally to decrease the line opacity, shifting the line Ty 
scale to the deeper hotter atmospheric layers and, thus, leading 
to weaker lines compared to the LTE case. 

3.2. NLTE effects in Sn\ 

The NLTE effects on the Srii levels are due to the non- 
equilibrium excitation i n the line transitions. Th i s was ex- 
tensively investigated bv pelvakova & Mashonkinal (Il997h and 
IShort & Hauschildtl(l2006h . and our results are qualitatively very 
similar to these studies. The analysis of radiative and collisional 
rates populating and de-populating the levels, as well as trial cal- 
culations with atomic models devoid of some key transitions, in 
particular the resonance and subordinate lines, allows us to draw 
the following conclusions. 

In the solar-metallicity models (Fig.|4] left panel), the ground 
state, 5s ^S (A^,, ~ 1130A) and the lowest odd level 5p^P° level 
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Table 3. NLTE abundance corrections for the 4607 Sr I line. Table 4. NLTE abundance corrections for the 4077 Sr II line. 



Tcff 
logg 


[Fe/H] 


2.20 


2.60 


3.00 


Anlte 
3.40 


3.80 


4.20 


4.60 


4400 


-3.0 


0.54 


0.51 


0.47 


— 


— 


— 


— 


4400 


-2.4 


0.49 


0.49 


0.49 


- 


- 


- 


- 


4400 


-1.2 


0.29 


0.29 


0.28 










4400 


-0.6 


0.27 


0.26 


0.24 










4400 


0.0 


0.19 


0.18 


0.16 


- 


- 


- 


- 


4800 


-3.0 


0.55 


0.55 


0.51 


0.49 


0.45 


— 


— 


4800 


-2.4 


0.50 


0.51 


0.51 


0.49 


0.46 


0.43 


0.36 


4800 


-1.2 


0.32 


0.33 


0.33 


0.33 


0.32 


0.30 


0.28 


4800 


-0.6 


0.17 


0.17 


0.18 


0.18 


0.17 


0.16 


0.14 


4800 


0.0 


0.19 


0.18 


0.18 


0.17 


0.15 


0.13 


0.10 


5200 


-3.0 


0.50 


0.52 


0.52 


0.50 


0.47 


0.45 


0.42 


5200 


-2.4 


0.47 


0.48 


0.49 


0.48 


0.45 


0.43 


0.41 


5200 


-1.2 


0.34 


0.35 


0.36 


0.36 


0.35 


0.34 


0.32 


5200 


-0.6 


0.19 


0.20 


0.20 


0.20 


0.20 


0.19 


0.17 


5200 


0.0 


0.09 


0.10 


0.10 


0.11 


0.11 


0.10 


0.09 


5600 


-3.0 


0.44 


0.45 


0.46 


0.46 


0.45 


0.42 


0.40 


5600 


-2.4 


0.40 


0.42 


0.43 


0.43 


0.43 


0.40 


0.38 


5600 


-1.2 


0.31 
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(/Ithr ~ 1540A , £■ ~ 3 eV) have nearly LTE populations out 
to log Tc ~ -3, and become underpopulated only in the outer- 
most layers, which are transparent to the radiation across their 
ionization edges. Still, there is small net radiative excitation in 
the transitions from these levels upwards, which very efficiently 
populates the levels with excitation energy E > 6 eV, such as 
5d^D, 6p^P°, and 7s -S (not shown in the plots) due to their 
small Boltzmann factors. The high-excited levels are, however, 
also over-ionized, which in turn, increases the number density 
of the Sriii ground state. The lowest metastable Srii level 4d-D 
(E ~ 1.8 eV) is nearly in thermal equilibrium with 5p^P° due 
to the dominance of collision rates in the deeper layers, however 
in the outer layers it gains appreciable overpopulation due to the 
photon losses in the near-IR transitions between the two levels. 
In fact, this process is a part of photon suction in a sequence of 
low-frequency transitions, which connect 4d to a large num- 
ber of high-levels close to the continuum. Many of these transi- 
tions become optically thin at logTc < -3, so that spontaneous 
de-excitation lead to the overpopulation 4d ^D, the lowest energy 
state of the cascading sequence. 

The departures from LTE change with decreasing metallic- 
ity (Fig. m right panel). First, the line transitions weaken and 
and over-population of 4d occurs in the deeper layers com- 
pared to the solar metallicity models. Also, radiative pumping in 
the resonance transitions at 4077 and 4215 A is more efficient 
due to larger radiative fluxes. This process leads to a marked 
over-population of the 5p^P° level at -1 < logTc < 0, espe- 
cially in the hotter model, Tefj - 6200 K. Outside these opti- 
cal depths, over-ionization dominates and all low-lying levels 
become underpopulated. The importance of strong line scatter- 
ing was tested by excluding the resonance lines and subordinate 
lines at 1 fim from the model atom. The overpopulation of the 



logg 


[Fe/H] 


2.20 


2.60 


3.00 


3.40 


3.80 


4.20 


4.60 


4400 


-3.9 


-0.07 


-0.04 


-0.04 


— 


— 


— 


— 


4400 


-3.0 


0.00 


0.00 


0.00 


- 


- 


- 


- 


4400 


-2.4 


0.00 


0.00 


0.00 


— 


— 


— 


— 


4400 


-1.2 


0.00 


0.00 


0.00 


- 


- 


- 


- 


4400 


-0.6 


0.00 


0.00 


0.00 










4400 


0.0 


0.00 


0.00 


0.00 


: 


: 






4800 


-3.9 


-0.11 


-0.08 


-0.04 


-0.01 


-0.01 


- 


- 


4800 


-3.0 


-0.02 


-0.02 


-0.02 


-0.01 


-0.01 


-0.01 


0.00 


4800 


-2.4 


0.00 


0.00 


0.00 


0.00 


-0.01 


-0.01 


0.00 


4800 


-1.2 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


4800 


-0.6 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


4800 


0.0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


5200 


-3.9 


-0.07 


-0.03 


0.01 


0.05 


0.08 


0.12 


0.13 


5200 


-3.0 


-0.05 


-0.06 


-0.06 


-0.05 


-0.04 


-0.03 


-0.01 


5200 


-2.4 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


5200 


-1.2 


0.00 


0.00 


0.00 


0.00 


-0.01 


-0.01 


-0.01 


5200 


-0.6 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


5200 


0.0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


5600 


-3.9 


0.05 


0.08 


0.11 


0.13 


0.14 


0.14 


0.15 


5600 


-3.0 


-0.07 


-0.08 


-0.09 


-0.08 


-0.06 


-0.04 


-0.01 


5600 


-2.4 


-0.01 


-0.03 


-0.03 


-0.04 


-0.04 


-0.04 


-0.03 


5600 


-1.2 


0.00 


0.00 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


5600 


-0.6 


0.00 


0.00 


0.00 


-0.01 


-0.01 


-0.01 


-0.01 


5600 


0.0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


6000 


-3.9 


0.19 


0.19 


0.19 


0.18 


0.17 


0.16 


0.15 


6000 


-3.0 


-0.04 


-0.06 


-0.06 


-0.05 


-0.04 


-0.02 


0.01 


6000 


-2.4 


-0.01 


-0.03 


-0.05 


-0.06 


-0.07 


-0.07 


-0.06 


6000 


-1.2 


0.01 


0.00 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


6000 


-0.6 


0.00 


0.00 


0.00 


-0.01 


-0.01 


-0.01 


-0.01 


6000 


0.0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



5p^P° and 4d^D vanished both for the metal-rich and metal- 
poor models. 

3.3. NLTE abundance corrections 

To quantify the effect of NLTE on the abundance determinations, 
we use the concept of the NLTE abundance correction Anlte, 
where: 

Anlte = log e(Sr)NLTE - log e(Sr)LTE ( 1 ) 

is the the logarithmic correction, which has to be applied to an 
LTE abundance determination of a specific line to obtain the 
exact value corresponding to the use of NLTE line formation. 
These values were computed for a grid of model atmospheres 
by equalizing the NLTE equivalent widths though varying the 
element abundance in the LTE calculations. The NLTE abun- 
dance corrections for the 4077, 4607 and 10036 line are given 
in Tables |3]|4] and |5] for a wide range of stellar parameters. Fig. 
ID (top panel) also shows the values of Anlte for selected model 
atmospheres as a function of effective temperature, gravity, and 
[Fe/H]. 

As evident from Fig. |5](top panel) and Table [3] the NLTE 
abundance corrections for the 4607 A are always positive, and 
they are maximum for cool giants with sub-solar metallicity, 
[Fe/H] < -1. This behavior can be easily understood by in- 
specting the plots of departure coefficients (Fig. H) top right 
panel): at the typical depths of the resonance line formation. 
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Table 5. NLTE abundance corrections for the 10036 Sr II line. 
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logTc ~ 0... - 0.5, the level populations are already depleted 
by a factor of two compared to LTE. In the models of warmer 
metal-poor dwarfs, the line formation is confined to the deeper 
layers, logTc ~ 0, where the NLTE effects in the Sri ioniza- 
tion balance are smaller. This is reflected in the behavior of the 
abundance corrections, which also decrease. Anlte are compara- 
tively small for the warmest models, corresponding to horizontal 
branch stars (Teff = 6200 and log^ - 2.2). Note, however, that 
at Teff > 6000 K and [Fe/H]< -1.5 the 4607 A line is very weak, 
Wx < 0.5 mA, and, thus, is not useful for abundance determina- 
tions. In contrast, the line could be suitable for the analysis in 
sufficiently good-quality (R > 20000, S/N >50) spectra of giants 

and sub-giants down to [Fe/H] 2.5. 

As described in the previous section, the formation of 
Srn lines is also significantly affected by the NLTE effects. 
Inspection of Table|4]and Fig.|5](middle panel) shows that for the 
resonance Sr ii line at 4077 A the NLTE abundance corrections 
are typically negative for stars with metallicity [Fe/H]> -3, but 
rapidly increase for more metal-poor stars with Teff > 5000 K 
and/or log g > 3. The amplitude of Anlte increases with Teff and 
decreasing log§. Another resonance Srn line at 4215 A shows 
the same pattern and the NLTE abundance corrections are nearly 
identical; we caution, however, that this line is blended by an Fe i 
line, and should be avoided in abundance analyses of cool stars. 
The implication is that, compared to NLTE, traditional LTE stud- 
ies relying on the resonance Sr ii lines under-estimate the abun- 
dance of Sr in very metal-poor stars, but somewhat over-estimate 



it for less metal-poor objects, such at that of thick and thin disk 
of the Galaxy. 

For the near-IR subordinate Srn lines, the NLTE abundance 
coiTections are typically negative (Table |5] and Fig. |5] bottom 
panel), being as large as -0.2 dex even for mildly metal-poor gi- 
ants; this behaviour is consistent with the discussion in Sect. l3.2l 
as the NLTE effects are due to strong line scattering. However, 
Anlte are usually within -0. 1 dex for dwarfs and become mildly 
positive only at very low metallicity, [Fe/H] < -3 and Teff > 
6000 K. The NLTE abundance corrections presented in Table 
|5] are rather similar for the other two members of the near-IR 
triplet, i.e., the lines at 10327 and 10914 A, and can be provided 
by request. 

Our NLTE abundance correcti ons are in line with the re- 
sults of IShort & Hauschildtl (l2006h . although they investigated 
the formation of the resonance Srn lines only. They find 
NLTE line strengthening for the giant models with [Fe/H] = 
-1 . . . - 2 and NLTE weakening for the models with [Fe/H] 

- - 4 ■ . ■ - 5. However, we note some qualitative differences 
with lAndrievskv et al.l (1201 ih . For example, for the near-IR Srii 
lines their abundance corrections at [Fe/H] - -3 and Teff = 
4800 - 5300 K are of the order of -0.3 to -0.5 dex and are very 
sensitive to surface gravity (their Fig. 7, top panel). In contrast, 
our values for the 10036 A line remain within ~ -0.1 dex for 
any log g value in this [Fe/H] and Teff range. For the resonance 
Sr II line at 4077 A, our Anlte are mildy negative for an y log g 
and Teff at [Fe/H] — -3, whereas fAndrievskv et al.l (1201 ih obtain 
large positive corrections for dwarfs and negative Anlte for gi- 
ants. The differences are most likely caused by the differences in 
the NLTE model atoms, as described in Sect. 12.21 

4. Application to the abundance analysis of 
late-type stars 

In addition to the Sun, three other late-type stars were selected 
for the abundance analysis. The principle selection criteria were 
that stellar parameters are well-constrained by independent tech- 
niques (interferometry, parallaxes), observations cover a large 
wavelength interval and their quality is very high for the 4607 
Sri line and the near-IR Srii triplet to be accurately measured. 
Stellar parameters for the warm turn-off star Procyon, moder- 
ately metal-poor giant Arcturus (HD 1 24897), and very metal- 
poor g iant H D 122563 were tak e n fromlKorn et"aD (j2003), Tsujl 
( 120091) . and iMashonkina et alj ( |2008|) . respectively (Table |6]l. 
The accuracy of these parameters was verified in our recent anal- 
ysis of Fe NLTE statistical equilibri um with mean 3D model at- 
mospheres dBergemann et al.ll2012l) . We stress that it is not our 
goal here to obtain Sr abundances for a large sample of stars, as 
this will be the focus of our forthcoming publications (Hansen 
et al. in prep., Ruchti et al. in prep.), but rather to test the perfor- 
mance of the new NLTE model atom. 

For the spectrum synthesis (Sect. 2.1), we used the so- 
lar KPNO flux spectrum (Kurucz et al. 1984). The spectra of 
Arcturus (HD 124897), Procyon (HD 61421), ar id HD 122563 
were retrieved from the UVESPOP database dBagnulo et al.l 
200l). These UVES (VLT, Paranal) spectra have a sUt- 
determined resolution of ~ 80000 and a signal-to-noise ratio 
SIN ~ 300 near 5000 A. 

The solar abundance was determined using five Sr lines 
from the Table [U adopting T^s = 5777, log^ = 4.44, [Fe/H] 

- 0, and = 1 km/s. The results are largely discrepant 
in LTE: logesri = 2.81 + 0.06 and logesrii = 3.06 + 0.1. 
The NLTE model atom successfully eliminates this problem. 
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Table 6. Basic parameters and the Sr abundances for the selected stars. 
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Fig. 6. NLTE and LTE line-by-Une Sr abundances for the se- 
lected stars. In contrast to LTE (red open circles), the Sri and 
Sr II lines provide consistent abundances in NLTE (black filled 
circles). 



The abundances from the two ionization stages are fully con- 
sistent logesri = 2.91 + 0.03 and logesiii = 2.94 ± 0.04 
in NLTE. Here, the error corresponds to one standard devia- 
tion of the line sample. The NLTE result is also in agreement 
with the meteoritic a bundance of Sr, loge = 2.90 + 0.03 dex 
jLodders et al.l |2009|). We note that the meteoritic abundance 
of Sr given bv lAsplund etaH (l2009h . Sr, log e = 2.88 + 0.03 
dex, is slightly lower, which is a consequence of different ref- 
erence Si abundances. The so l ar Sr a bundance was also deter - 
mined by iBarklem & O'Maral (l2000h . lAndrievskv et al.l (1201 ih . 
and Mashonkina & Gehren (2001). The last two references per- 
form NLTE calculations, which are generally in agreement with 
our results, both in terms of NLTE abundance corrections and 
absolute abundances. We caution, however, that this agreement 
is not particularly telling, because other parameters in the model- 
ing are different. Mashonkina & Gehren (2001) adopted = 0.8 



km/s and employed the older (ODF) ver s ion of the MAFAGS 
model atmospheres. iBarklem & O'Maral ( |2000|) recover mete- 
oritic Sr abundance from the neutral Sr lines even in LTE, which 
is not unexpected as they used = 0.85 km/s and the solar 
Holweger-Miiller (1974) model atmospher^ 

The LTE and NLTE Sr abundances determined for the refer- 
ence stars are given in Table |6] and a line-by-line comparison 
is shown in Fig. |6] Few examples of synthetic and observed 
line profiles are shown in Fig. |7] Comparison of the [Sri/Fe] 
and [Srii/Fe] from Table |6] reveals that the LTE approach pre- 
dicts a systematic discrepancy between the lines of two ioniza- 
tion stages, which is also evident from the Fig.|6] The difference 
between the Sr i and the near-IR Sr ii lines, which are almost in- 
sensitive to damping, is of the order 0.3 to 0.5 dex. This is far 
beyond the uncertainties in stellar parameters, which are ~ 80 
K for Jeff, 0.1 dex for \ogg, and 0.1 dex for [Fe/H]. Also, the 
uncertainties of the transition probabilities are far smaller than 
that (Sect.l23Tl. 

The NLTE model atom provides, in contrast, realistic ion- 
ization balance, significantly reducing the scatter between the Sr 
lines of different ionization stages and excitation energies (Fig. 
|6|l. In this respect, especially important is a good agreement be- 
tween the NLTE abundances from the near-IR Sr ii and visual 
resonance Sri line, which are very sensitive to NLTE and the 
effects are of different nature, i.e., overionization-stipulated de- 
crease of opacity in the 4607 line and source function depletion 
in the lines of the lyum triplet. 

The analysis of a larger sample of stars supports the impor- 
tance of NLTE effects in the ionization balance of Sr. This is ev- 
ident from Fig. [8] which shows the difference between the mean 
abundances computed using the Sri and Srn lines for a sample 
of thick-disk and halo RAVE stars with high-resolution spectro- 
scopic follow-up observations (see Ruchti et al. 2011). In LTE, 
the mean offset between the two ionization stages is ~ 0.4 dex 
and it is alleviated when the NLTE effects are taken into account. 
A detailed analysis of the data will be presented in Ruchti et al. 
(in prep.). 

We thus conclude that NLTE must be taken into account in 
abundance analysis of resonance Sr i and subordinate Sr ii lines 
in spectra of late-type stars. Also, at low metallicities, the reso- 
nance Srii lines are sensitive to NLTE effects, and abundances 
determined in LTE are largely under-estimated. It is also impor- 
tant to note that the NLTE effects for Sr might be larger than our 
estimates, because of interlocki ng with NLTE-affec t ed Un es of 
other atoms, as demonstrated bv lShort & Hauschildll (l2006l) . 

5. Conclusions 

We investigate statistical equilibrium of Sr in the atmospheres 
of late-type stars. The NLTE model atom was constructed us- 
ing new atomic data, computed with the R-matrix method. 



' The semi-empirical Holweger-Miiller model atmosphere is hotter 
than all solar theoretical models by ~ 150 K over the line formation 
layers and, in this way, it 'mimics' NLTE effects for some atoms 
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Fig. 8. NLTE and LTE abundance differences between Sri and 
Srii lines for the representative sample of thick-disk and halo 
stars from Ruchti et al. (201 1). 



which include levels, transition probabilities, photoionization 
and electron-impact excitation cross-sections. 

The NLTE effects are significant for the Sr i resonance line 
at 4607 A and are stipulated by overionization. For the reso- 
nance Srii lines, typically observed in metal-poor FGK stars, 
our model predicts small deviations from LTE, although NLTE 
effects become important at very low metallicity, [Fe/H] < -3, 
reaching +0.2 dex for very metal-poor turnoff stars. The near-IR 
Sr II triplet shows substantial sensitivity to NLTE, and the NLTE 
corrections vary between -i-O.l and -0.2 dex depending on the 
model metallicity, temperature, and surface gravity. 

The NLTE model atom was applied to the analysis of the 
Sun and three stars with well-constrained stellar parameters. In 
contrast to LTE approach, the NLTE model recovers ionization 
and excitation equilibria of Sr for all these reference stars, thus 
confirming that the NLTE modeling approach developed in this 
work provides a solid base for future abundance determinations 
of late-type stars. 

The NLTE abundance corrections are provided for the im- 
portant diagnostic Sri and Srii lines for a grid of model at- 
mospheres in the following range of stellar parameters 4400 < 
Tetf < 6000 K, 2.2 <logg < 4.6, -4 < [Fe/H] < 0; for other 
Sr lines and other models the abundance corrections can be also 
computed by request. 
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Fig. 5. NLTE abundance corrections for the Sri and Srii lines. Note different x-axis scales. 
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